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Preface 
This  report  contains  two  parts.   Part  1  derives  the 
normal  mode  stability  equation  for  a  diffuse  pinch  of 
guiding  center  plasma  and  examines  the  guiding  center 
Suydam  instability  and  its  growth  rates.   Part  2  studies 
the  pinch  stability  problem  in  the  low  shear  limit. 
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Part  1 
STABILITY  EQUATION  AND  LOCAL  INSTABILITY 
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I .    Introduction 

A  convenient  approach  in  studying  plasma  stability 
problems  is  a   6W  variational  analysis.   This  yields  a 
6W  stability  criteria  but  gives  no  information  on  the  growth 
rate.   An  alternative  approach  is  to  use  the  normal  mode 
analysis  for  the  full  set  of  equations  of  motion,  or,  equiva- 
lently  the  variational  analysis  including  the  kinetic  energy. 
This  not  only  gives  the  stability  criteria  but  also  yields 
information  regarding  the  growth  rate  and  the  spectrum  of 
the  problem. 

The  present  work  is  concerned  with  the  stability  problem 
of  a  cylindrical  linear  diffuse  pinch  of  collisionless  guiding 
center  plasma.   For  this  problem,  the  6W  variational  analysis 
has  been  studied  by  Aleksin  and  Yashin.    The  purpose  of  the 
present  work  is  to  study  this  problem  by  using  the  normal 

mode  analysis  based  on  the  full  set  of  equations  of  motion  of 

2  3  4 
the  lowest  order  guiding  center  theory  '  '   and  obtain  informa- 
tion regarding  the  growth  rate  and  the  spectrum  of  the  problem. 
It  is  found  that  the  full  set  of  equations  of  motion  can  be 
reduced  to  a  single  second  order  ordinary  differential  equa- 
tion [see  Eq.  (57)]  which  is  convenient  for  both  analytic  and 
numerical  purposes.   This  second  order  stability  equation, 
which  contains  the  growth  rate  w  in  it  and  which  reduces  to 
the  equation  of  Aleksin  and  Yashin  in  the  special  limit  w  =  0, 
can  be  viewed  as  the  guiding  center  version  of  the  Hain-Ltist 
equation  for  the  MHD  diffuse  pinch^.  It  is  somewhat  similar 


to  the  Hain-Ltist  equation  except  that  its  coefficients  contain 
Cauchy-type  integrals  of  the  equilibrium  distribution  functions. 

In  a  first  application,  this  equation  is  then  used  to 
study  the  Suydam-type  instability  and  the  associated  spectrum. 
The  result  confirms  the  MHD  result  of  Grad  that  Suydam  in- 
stability has  infinitely  many  discrete  unstable  eigenvalues 
with  the  origin  as  accumulation  point.   And  it  is  shown  that 
the  eigenvalues  w   for  large  n  satisfy  the  linear  relationship 

u  .t/o)   =  exp(-Tr/2c) 
n+i   n 

where  c  is  a  constant  parameter  [see  eq.  (74)];  this  result 
has  been  born  out  by  the  numerical  calculation  of  Goedbloed  and 
Sakanaka   for  the  MHD  diffuse  pinch. 

Further  stability  study,  including  line-preserving  and 
semi-local  instabilities  ,  based  on  the  stability  equation 
[eq.  (57)]  is  currently  in  progress. 

It  is  assumed  that  in  the  equilibrium  state  the  plasma 
properties  and  the  magnetic  field  vary  with  the  radial  variable 
r,  but  are  independent  of  the  azimuthal  variable  (J)  and  the 
axial  variable  z.   The  equations  of  motion  are  studied  in  Section 
II  and  it  is  shown  in  Section  III  that  the  result  can  be  reduced 
to  a  second  order  ordinary  differential  equation.   The  boundary 
conditions  are  formulated  in  Section  IV  by  assuming  that  the 
plasma  is  either  bounded  by  a  conducting  wall  or  is  separated 
from  the  wall  by  a  vacuum  region.   In  Section  V,  the  eigenvalue 
problem  associated  with  the  Suydam-type  instability  is  studied 
and  the  spectrixm  is  also  discussed. 


II   Basic  Equations 

Consider  a  plasma  consisting  of  electrons  and  one  kind 

+ 
of  positive  ion  with  distribution  functions  F  .   We  shall  use 

the  version  of  guiding  center  theory  given  in  Reference  2  and 

write  F   =  F  ((-l,v,x,  t)  where  \i   is  the  magnetic  moment  and  v  is 

the  particle  velocity  along  the  magnetic  field  line.   In  terms 

±  +  ±  ± 

of  P  ,  the  density  p—  and  the  pressure  components  P.,  and  Pp 


+         r         2  ± 

[p,P^,P2]-  =  ^[l^v  ,nB]F  dudv  (i) 


and  for  the  plasma  as  a  whole  we  have  p  =  p  +  p~,  P  =  P  +  p' 
and  Pp  =  Pp  +Pp.   The  pressure  tensor  is  then 


(2) 


o 
where  A  =  (P  -  Pp)/B  and  I  is  the  unit  tensor. 

Let  U  "be  the  macroscopic  plasma  velocity  and  introduce 

Ji  "^  y  "  (JJ'6)S  where  g,  is  the  unit  vector  g/B.   The  basic 

equations  are 

dU 
P  di  ="   -^^+(^^g)  ^g  '  (3) 

SB 

^=  Vx(uxB)  ,  (4) 


V-B  =  0  .  (5) 

On  the  other  hand  the  motion  of  guiding  centers  along  the 


magnetic  field  lines  is  governed  by  the  one-dimensional  kinetic 
equation 


0  =  fl   +  V.[(u  +vg)F  ]  +  K.u(vF  ''v  +  [a  T  +£.V(Au^  -M.B)]F  ^   (6) 


where  the  subscript  v  denotes  differentiation  with  respect  to 
V,  K  is  the  magnetic  field  curvature, 

T  =  P'[(VP/p)  +  -  (VP/p)-]  (7) 

+ 
and  the  constants  a   is  given  by 

a-  =  (e/M)7f(e/M)+-  (e/M)"]  ,  (8) 

+      + 
where  e   and  M  are  respectively  the  charge  (signed)  and  the 

mass  of  the  ions  and  electrons. 

Following  the  nonnal-mode  method,  we  consider  the 

perturbations 

F  =  F°  +f  exp(ie)  ,  (9) 

^1,2  =  Pl,2+Pl,2  e^P(iS)  '  (10) 

B  =  B°  +b  exp(i9),  (11) 

U  =  U  exp(i9),   p  =3  p°+p^  exp(i0)  ,  (12) 
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where  9   =    (-cot  +  kz  +m<t))  in  the  cylindrical  coordinates  (r,<l),z) 
and  the  superscript  o  denotes  equilibrium  values.   Only- 
unstable  modes  are  considered  and  it  is  therefore  assumed 
throughout  that  w  is  not  real.   For  the  equilibrium  state,  the 
pressure  balance  requires  that 


P.+^"^(°B°)2  =   0  (13) 


jhere   P^  =   Pp  +^(B°)    ,    the   dot   denotes   differentiation  with 


respect   to   r,    o    is   given  by 


r^   =  1-A°  (14) 


and  A°  is  defined  in  (2) .   It  is  further  assumed  that  for  the 
equilibrium  state  both  o      and  the  quantity  K-|^(0)  given  by  (38) 


'^    >  0,      K^(0)  >  0  (15) 


so  that  the  fire  hose  and  mirror  instabilities  "-^  ■'  associated  with 
a   <  0  and  K  (0)  <  0  will  not  enter  the  present  pinch  problem. 

For  convenience  of  writing,  we  shall  introduce  the 
following  notation 


g  =  B°.b,   b^  =  (B°b^-B°bJ/B°  (16) 


a  =  div  u,   u^  =  (B°u^  -B°uJ/B°  (1?) 
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T  =  P2  +g.   C  =  P^  -P2  (18) 

H  =  kB^+(m/r)Bj,   v  =  B7(rB°)  ,  (19) 

E  =  kB°-(m/r)B°  ,  (20) 

and  denote  by  r  and  $  the  unit  vectors  in  the  r  and  * 
directions  respectively. 

Now  substituting  (10) -(12)  into  the  equation  of  motion 
(3)  yield  the  linearized  equation 

-p  iCDU  =  -Vt-b"^(C  -2Ag)(iHB-r"^B^?)  +h^[  0^(6  +  r"^B^$)  -  AB] 

(21) 
+  a2[iHb-2r-\b^r] 

where  we  have  omitted  the  superscript  o  from  the  equilibrium 
values  p°,  B°  and  A°  for  simplicity.   Similarly,  equations  (4) 
and  (5)  yield 


icDb  =  -  Hu  +ag+u  (B- r"^B^$)  ,  (22) 


b^+(b/r)  +ikb^  +i(m/r)b^  =  0  ,  (23) 

where  (23)  is  not  independent  of  (22). 

Now  substituting  (9) -(12)  into  (6)  gives  the  following 
linearized  kinetic  equation 


+         _T  +   +     _2   ±  -P   ±        1  .+     T  + 

f   =  IH   B-'F^Lari-B   (a  V  +  B^l)g  ]  &  +B   gF  +(la))^(F  -  B"-^BF  )u 


-  u^F^[(rB2)"^B^v+H(Hi3)"^(a  S  -  nB)  ]  6 


where  again  we  have  omitted  the  superscript  o  from  F  and  B. 
The  function  r],  5,  V  and  S  are 


r 
(24) 


n  =  (p/p)^-  (Pi/p)"  ,  (25) 


=  (icD-  lHB"^v)"^  (26) 


V=  [{F^-F^)/pr-[{P^-F^)/p]-    ,  (27) 

S  =  (P/p)+-  (P/p)--  V(B/B)  ,  (28) 

Where  it  is  recalled  that  p  is  the  equilibrium  plasma 

density,  and  it  is  noted  that  5  Is  always  finite  since  only 

+ 
non-real  values  of  w  are  considered.   The  perturbation  p.,  and 

+  f- 

P2  are  related  to  the  perturbed  distribution  function  f  by 

±   r  2  i 

p^  =  /  V  f  dudv  ,  (29) 

P2  =  Bj  nf   dUdv+B  ^P2g  ,  (30) 

Pi  =  P^  +Pl  ,   P2  =  P2  +P2  •  (31) 


The  basic  equations  are  eqns.  (21)-(24)  while  eqns.  (I6)- 
(20)  and  eqns.  (25)-(31)  furnish  the  definitions  of  the  various 
quantities  that  are  involved  in  the  basic  equations. 
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Ill  The  Stability  Equation 

It  Is  possible  to  reduce  the  basic  equations  (21) -(24) 

in  the  last  section  to  a  single  second  order  ordinary 

differential  equation.   This  will  be  done  in  the  present 

section.   The  derivation  involves  considerable  arithmetic 

manipulations  and  a  complete  description  of  the  derivation 

will  not  be  presented.   In  what  follows  only  the  major  steps 

in  the  derivation  will  be  given. 

+ 
First  substitute  expression  (24)  for  f  into  (29)  and 

(30)  and  use  (25)  and  (31)-   This  then  gives  rise  to  three 

equations  connecting  T],V-i>Vn)  S>   ^-i^d  u  .   We  can  then  solve 

these  equations  for  y],  p,  and  p^  in  terms  of  g  and  u  .   This 

result  then  enables  us  to  express  p,  and  Pp  in  terms  of 

linear  combinations  of  g  and  u  .   If  these  expressions  are 

used  in  (l8)  for  x  and  C,  we  can  then  write 

T  =  K^g  +  {lG,)-\n^   ,  (32) 

C  =  L^g  +  (icD)"%u^  .  (33) 

The  coefficients  K-,,Kp,L  and  Lp  can  be  expressed  in  a 
convenient  form  if  we  introduce  the  functions 

!<*'")  =/"v^/"-*  [(^)"  +  (-l)"{^)-l  d.  .  (34) 

"^co      *-^  o       p  p 

J^(w)  =  1(2,0)  -  [I(1,1)I(1,1)/I(0,2)]  (35) 

-'■■       ,9.2^^)  =   l(l'O)  -[I(0,l)l(l,l)/I(0,2)]  (36) 
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^^(vj)    =  1(0,0)  -  [1(0,1)1(0, 1)/I(0,2)]  (37) 

w  =  ooB/H  (38) 

In  terms  of -JJ,,  ^p  and  ^fi^,  the  coefficients  in  eqns .  (32) 
and  (33)  can  be  written  as 

K-^(w)  =  2B"^P2  +l+-^^(w)  ,  (39) 

K2(w)  =  -BBK^(w)  -  r"^B^[a2  _  b~\^--^^{v^)  ]  ,  (40) 

L-^(w)  =  a2+2A-K^(w)  +W^B"^.'^2(^)  '  (^1) 

L2(w)  =  P-L  -BBA-K2(w)  -w^B^)^2(^)  -  ^~"^(B'*'/B)^[pw^- ^/^^^(w)].  (42) 

The  main  result  from  the  one-dimensional  guiding  center 
kinetic  equation  is  given  by  the  expressions  (32)  and  (33), 
which  is  all  that  is  needed  in  the  equations  of  motion.   In 
this  sense,  the  guiding  center  equation  along  the  field  plays 
the  role  of  constitutive  equations  which  furnish  the  relation 
between  t, C  and  the  basic  variables  g  and  u  . 

We  now  turn  to  the  equations  of  motion.   The  r-component 
of  (21)  is  simply 


-pia)u^  =  -T+r-l(B/B)2(C-2Ag)  +a2tiHb^-2r-\b^^  ^       ^^^^ 


In  the  direction  perpendicular  to  both  r  and  B,  U  and 
u  have  the  same  component  u  .   The  equation  for  u  can  be 
derived  from  (21)  and  the  result  is 


■  11- 


•plcoBu^  =  iET+ia2Hb^+b^o2B^B^[(C/v)  +(2/r)]  .         (44) 


where  v  Is  given  in  (19). 

On  the  other  hand,  the  divergence  of  u  can  be  written  in 
terms  of  u  and  u  as 


a  =  u^  +  r"^u^  -  iB"^E  u^  (45) 


and  from  the  magnetic  field  equations  (22)  and  (23),  it 
follows  that 


iEBu^  =  -ia)g  +  ( BB  +r"-'-B^)u^  +  B^u^  ,  (46) 


l^b^  =  B^B^(v/v)u^-iHBu^  ,  (47) 


i^b^  =  aB^  +  (B^-r-lB^)u^-iHB^B-\  ,  (48) 


b^  =  -m-\   ,  (49) 


from  which  it  is  readily  seen  that  u  , b  ,b.  and  b   can  all  be 
expressed  in  terms  of  u  ,u  and  g. 
into  (44)  and  using  (32),  we  obtain 


expressed  in  terms  of  u  ,u  and  g.   Substituting  these  expressions 


a.cDg  =  NB^u^  +Q  u^  (50) 


where  N,0  and  Q  are  given  by 


N  =  fyu^  -  a^H^  ,  (51) 


■12- 


N(BB+r"^B^)  +^K^-2t~^    o-b^B^EH  ,  (52) 


O  --  N  -  S^K-L  .  (  53) 


Finally,  we  substitute  (32),  (33),  (48)  and  (49)  Into 
(43),  using  (46)  and  (50).   The  result  is  a  second  order 
ordinary  differential  equation 


d  ,^1 


K,NB^  da  du 


[^O-  -d?]  +  ^l(^'  ^)  -^  +  ^o^^'  r)  u^  =  0         (  54) 


(55) 


where  the  coefficients  C   and  C,  are 
o      1 


C-^(w,r)  =  0"^(K^Q+K20)  -  (rO)"^(L^-  2A)NB^ 


+  2(rO)  ^  a  kEB  K-^+2r    a^B^ 


C^(w,r)  =  N +-^[0"^(K^Q+0K2)  ]  +2r"-^o2] 


2(rE)"^a2kB^[(Q/0)  -BB- r"S^]  (56) 


r-l(B/B)2  [L2+(Q/0)  (a^-K^)] 


A  further  simplification  for  (54)  is  possible.   For  this 
purpose,  we  substitute  (39)-(4l)  into  (55)  and  after  some 
manipulations  it  is  found  that 


C^  =  K-|_NB^(rO)"-'- 


■13- 


and  as  a  consequence,  eqn.  (54)  may  be  written  as 

-2-[ — k ^]  +C  (w,r)ru  =  0  ,  (57) 

dr"-   O     dr-^   o^    '   r  ^^'  ^ 


where  C   is  given  in  (se).   Eqn.  (57)  Is  the  stability 
equation  which  forms  the  basis  for  all  subsequent  discussions. 
In  the  limiting  case  for  w  =  0,  eqn.  (57)  agrees  with  the 
variational  result  obtained  by  Aleksin  and  Yashin  . 

It  is  also  noted  that  if  F—  are  symmetric  in  v  and  non- 

2 
increasing  in  v  ,  then  K, (0)  >  0  implies  that  K  (W)  >  0  for 

all  imaginary  W.   It  then  follows  that  (57^  is  non-singular. 


ih- 


IV   Boundary  Conditions 

If  the  plasma  fills  the  cylinder  up  to  the  conducting 
wall  at  r  =  R  ,  then  the  boundary  condition  for  (56)  at 
r  =  R  is  simply 


u^  (R^)  =  0  .  (58) 


On  the  other  hand  if  the  plasma  in  equilibrium  has  a 
vacuim-plasma  interface  at  r  =  R  ,  then  the  interface  must 


also  be  perturbed  as 


r  =  R^  +  4  exp(i9) 


and  the  kinematic  boundary  condition  at  the  interface  is 


103^  +u^  =  0  .  (59) 


Outside  the  plasma,  the  vacuum  magnetic  field  may  be 
denoted  by 


B  =  B°+b^  (60) 


and  the  perturbed  external  field  b   can  be  expressed  as 


bg  =  V)^  (61) 


-15- 


where  x  =  x(^)  exp(ie)  and  x(^)  is  a  linear  combination  of 
the  modified  Bessel  functions  of  order  m. 

At  the  interface,  continuity  of  the  normal  component 
of  the  magnetic  field  requires 


b^-i(Hg-H)^  =  -X(r)  ,  (62) 


where  H  =  kB   +(m/r)B^  ,  and  we  have  dropped  the  superscript 
o  from  B°.   The  pressure  balance  gives 


^[-  ^'B^.B,3B,^]-fT  =   H^X.  (63) 


If  we  substitute  (32)  for  t,  (48)  for  b^  and  (59)  for 
i   into  (62)  and  (63),  we  find  that 


-oox  =  (Hg-2H)u^  ,  (64) 


b2 

-CDH^X  =  (K3  -  -|^)u^  +0-\nb2u^  ,  (65) 


where  K^  is 


K^  =  {K^(^)/a+K^  +  a%^^    .  (66) 


Equation  (64)  and  (65)  can  be  rearranged  to  give  the 
final  boundary  condition 


-16- 


b2 
0  =  O-^K^NB^  u^+[K3--^  +  (2H-H^)H^(x/x)]u^  (6?) 


at  r  =  R  ,  where  x/x  is  completely  determined  by  the 
position  of  the  conducting  wall  at  r  =  R,  which  bounds  the 
vacuum  region. 

Finally,  it  Is  noted  that  r  =  0  Is  a  regular  singular 
point  of  fS?).   It  is  not  difficult  to  show  by  displaying 
the  indicial  equation  that  as  r  -*  0,  one  solution  of  (57j  tends 
to  zero  while  any  other  linear  independent  solution  diverges 
in  such  a  way  that 


/ 


(u^)  r  dr 


is  infinite  which  is  unacceptable.   Therefore  the  boundary 
condition  at  r  =  0  may  be  taken  as 


u^(0)  =  0.  (mj^O)  (68) 


-17- 


V    Suydam  Instability  and  the  Associated  Eigenvalue  Problem 

In  this  section  we  consider  the  case  In  which  H  vanishes 

at  a  point  r  =  r  In  the  plasma,  and  we  shall  be  concerned 

with  small  values  of  w.   The  eigenvalue  problem  to  be 

considered  consists  of  the  second  oi^der  ODE  In  (57)  and  the 

boundary  conditions  which  may  be  e^'.ther  (67)  and  (68),  or 

(58)  and  (68). 

To  facilitate  our  discussions,  it  is  expedient  to 

consider  three  Intervals  (o,  r  -  e),  (^q  "  e-ij  ^q  +s-|)  sind 

fr  +  e,  R  )  where  £<<£-,  <<  r  ,   In  the  two  side  Intervals 
^  o      o'  1     o 

(o,r  -  e)  and  (r  +  e,  R  ),  H  is  bounded  away  from  zero  so 
that  pw   is  negligible  compared  to  o  IT  and  eqn.  (57)  may  be 
simplified  as 


2..2„2 


du 


,   rK,  0  cr^H  B' 

^[-^4 2 d?]  +rC  (0,r)u  =  0  .  (69) 

^^  o^Yf   +E^K^(0)   dr     o      r 

As  r  in  the  two  side  intervals  (o,  r  -  e)  and  (^  +  ^>^q) 

tends  to  r  +  e,  we  have 
o  — 


H   =  m  B  V  X  (70) 


-k/m  (71) 


where  x  =  ( r -r  )  and  after  some  computations  we  find  that 


2  2 

^o^o^^o)  =  "T^t^ +Po+^(i-Trw)J  (72) 

o 


5j[P,.P2.^(l-iJ^ 


Using  the  variable  x  =  r  -  r  ,  we  can  write  eqn.  (69)  as 
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^4^-^"r  =  0  (73) 

2 
where  the  coefficient  q  is  C  (0,r  )  divided  by  (chb/E)  . 

The  solution  of  this  equation  is  oscillatory  if  q  >  jr,    given 

by 


xl-l/^  3in  @.   1x1-1/2  ^^^ 


2 
where  0  =  c  log  x  and  c  is  the  index 


=  =  i(5-i)^/2.  "^' 


It  may  be  claimed  that  a  sufficient  condition  for 
instability  is  that  q  >  jr,  which  may  be  written  in  full  as 


i/^)'^^^  -   K^  ]  "-    ^^^b2(v/v)2  <  0  .        (75) 


This  is  the  guiding  center  version  of  the  Suydam 
condition  for  instability.   It  may  be  noted  that  if  the 
electrons  and  the  ions  have  the  same  distribution  function 
and  if  the  distribution  function  is  a  two-temperature 

Maxwellian,  then  K-,  (0)  =  1.   If,  moreover,  the  pressure  is 

2 
isotropic,  then  o   =  i.   Under  these  circumstances,  (75) 

9 

reduces  to  Suydam 's  result  . 

Condition  (75)  has  previously  been  obtained  by  Aleksin 
and  Yashin-'-  using  a  variational  analysis. 

To  show  that  (75)  indeed  implies  instability  and  to 
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examine  the  nature  of  the  eigenvalue  problem,  we  now  turn  to 

the  center  layer  r  -e,  <r<r  +£.,. 
"^     o   1  —   —  o   1 

We  only  need,  to  deal  with  purely  imaginary  od  and  can 
therefore  write 


p(r^)a)2E-2  =  -A^  .  (76) 


In  the  center  layer,  we  return  to  the  original  equation, 
( 57) ,  and  introduce  a  new  variable 


'x  =  x/A 


in  terms  of  which  eqn.  ( 57)  may  be  written  in  the  form 


o  dli  5 

[rC^  )-^]  ^^r.^X^'^J'^^  =   0  (77) 


d^'-  ^   '  -.t^^        00^   '  o'  r 


valid  in  the  center  layer  r  -e-,_<r_<r  +e. 


r  +  e^  .   The  function 


rCx^)  is 


rCx^)  =  r  (B/E)^[l  +m^B^v^^] 


where  B,  E,  B  and  v  are  evaluated  at  r  =  r  . 
'    '    z  o 

The  advantage  of  using  x  is  that  eqn.  (77)  is  now 
independent  of  w.   As  x ->•  +  00 ,  eqn.  (77)  has  the  same  form 
as  eqn.  (73)  and  therefore  the  solutions  of  (77)  have  the 
asymptotic  behavior 
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l-^r^^^  sin  -O;   1-^1-1/2  cos  '&  (78) 


as  "x  -*  +  oo ,  where  ®  =  c  log  'x  . 

Since  the  coeffifients  of  (77)  are  symmetric  in  '^, 

we  can  speak  of  a  symmetric  solution  u  and  an  ant i- symmetric 

solution  u  .   According  to  (78)^  we  may  express  the  asymptotic 

behavior  of  u  and  u^  as 
S       a 

u^   =    \3i\~^^^   sin(@  +  |),     'x-^  +  oo,  (79) 

vi^  =  ±\'ir-^^^[oos{&  +  ^)  +A   sln{^  +1)],      5c  ^  +  00  (80) 

where  the  constants  1  and  A  are  independent  of  o). 

Now  it  is  possible  to  put  the  pieces  together.   First 

we  start  with  a  solution  u  of  (57)  which  satisfies  the 

boundary  condition  (58)  or  (67)  at  r  =  R^.   As  r  moves  from 

R  to  r  +  e,  this  solution  becomes 
o     o   ' 


u^  ^  \x\-^^^   sin(0  +$^)  .  (81) 


And  on  the  other  hand,  in  the  interval  (0,r  -  e),  we  start 
with  the  boundary  condition  (68),  and  as  r -*  r^  -  e,  we  have 


u^  ^  gIxI"^^^  sin(@  +1_)  .  (82) 


It  is  noted  that  I  are  also  independent  of  03. 

What  remains  to  be  shown  is  that  it  is  always  possible 
to  find  a  real  A  and  a  solution  u  of  (77)  such  that  u^  matches 
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(8l)  as  ^-*  +00  and  (82)  as  ^  -*  -od  . 

To  illustrate  the  basic  ideas  in  the  proof,  let  us  first 
consider  the  special  case  in  which  I"  =  l_(mod.  2Tr) .  In  this 
case,  we  only  need  to  take 


^  //^^N     T,-l/2 


which  goes  to 


Tl^ -*  \'x\\~^^^   sin(0-c  log  A^+f) 


as  "S! -*  +  00  .   In  order  for  this  solution  to  match  (81)  and  (82) 
it  is  only  necessary  to  choose  A  such  that 


c  log  A  =  i  -  i  -  7m 


where  n  is  any  integer.   We  therefore  have 


A^  =  -p(r^)^/E-2 


(83) 


=  exp[c~^(i-  i^-7m)  ]  . 


If  i  =  1_  +ir,  the  same  argument  goes  through  by  taking  G  =  -1. 

It  is  seen  from  (83)  that  if  n  is  large,  then  A  will  be 
small  which  justifies  our  consideration  that  w  is  small. 

Now  we  proceed  to  the  general  case  in  which  J  /   $_  (mod.  tt) 
We  start  with  (81)  and  in  terms  of  the  variable  "4,    eqn.  (81) 
becomes      '     •'■'-■''■ 
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u^  =    |ScAr^^^[cos   ^^  sln(^+I)  +sln  ^^  cos(e+i)]    ,  (84) 


where   the   constant   ip     is   given  by 


^^  =  c   log  A^  +5^-1   .  (85) 


The  expression  in  (84)  can  be  decomposed  into  a  symmetric 
part  given  by  (79)  and  an  anti-symmetric  part  by  (80)  and  the 
result  is 

u^  =  |5?'Ar-'-^^{(cos  ^^-A  sin  ^_^)sin('©  +1) 

(86) 
+  sin  ^_^[cos(e  +i)  +A  sin(^  +1)  ]]  . 


Now  in  order  for  u  to  match  (86)  as  "x -►  +00,  we  choose 


^     -1/2 

u^=x   ^  [(cos  ^^-A  sin  ^^)Ug  +(sin  ^^)u^]  .         (8?) 


It  then  follows  from  (79)  and  (80)  that  as  'x -^  -00 


u^  =  |xAr^^^[(cos  ^^-2A  sin  ^_^)sin('@  +£) 


-(sin  Tp   )    cosC^H-^)]  . 

which  must  match  (82).   To  see  that  this  is  possible,  we 
rewrite  (82)  in  terms  of  0  as 

u^  =  G|xAr-'-'^^[cos  ^  sin(^ +  i)  +sin  ^  cos(^+l)]        (89) 
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where  ^  can  be  found  in  (85).   Now  in  order  for  (88)  and  (89) 
to  match,  it  is  sufficient  to  have 

sin  i//_(cos  ^_^-2A  sin  ^_^)  =  -sin  ^_^  cos  Tp_ 

or,  upon  using  (85), 

sin  r-A[cos  r-cos(I^-I_)]  =  0  (90) 

where  r  stands  for 


r  =  Tp_^+Ti/_ 


(91) 


2c  log  A^  +1_   +1^-  21. 


The  following  argument  shows  that  (90)?  a-s  an  equation 
for  r,  has  infinitely  many  solutions.   Let  us  assume  that 
cos(I  -^  )  >  0.   We  choose  T-.    and  T^^   such  that 


cos  r 


3. 


^^2  =  cos(i^-  i_) 


and  0  <  r^  ^  J  and  -^Tr  ^  r^    <   2Tr.   Then,  for  r  =  r^,  sin  r 
is  positive  and  therefore  the  left  hand  side  of  (90)  is 
positive,  and  for  r  =  F^,    it  is  negative.   Therefore  eqn. 
(90)  must  have  a  root  r^  in  (r^,r2).   If  cos(i^-i_)  <  0,  we 
only  have  to  consider  X  _<  r-,  <  fr  and  tt  <  r^  j^  ■2'"'- 
The  general  solution  of  (90)  is  then 
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r^  +  2n7r 


for  n  being  any  integer,  and  consequently 


a2  =  .p(r^)a.V2 


(92) 


exp[c"^(5-|l^--|$_  +-|r^-n7r) 


2^+  2*-   2' 

and  only  large  n  will  be  taken  so  that  cd  is  indeed  small. 

From  (92)  and  (83),  it  is  seen  that  there  are 
infinitely  many  eigenvalues  with  cd  =  0  as  an  accumulation 
point.   This  is  in  agreement  with  the  magneto-hydrodynamic 
result  obtained  by  Grad^.     On  the  other  hand,  if  c  is 
small,  the  maximum  growth  rate  is  bounded  by 


|a)|^<  [E2/p(r^)]exp(-iyc) 


where  1^  is  obtained  as  the  smallest  possible  positive  value 
of  nTT +-^i  +  2-.f_  -  Q-r^  -  I.  for  all  possible  integer  n. 

It  is  also  noted  from  (87)  that  the  amplitude  of  u   is 

-1/2 
by  a  factor  A     larger  than  the  amplitude  of  the  eigenfunction 

in  the  side  intervals  (0,r  -  e)    and  (r  + e,  R  ) .   This  indicates 
^  '  o    ^      ^  o     o^ 

that  the  eigenfunction  is  localized  in  the  center  layer  around 
r  =  r  ,  while  the  thickness  of  the  layer  is  of  the  order  Ar  . 
Finally,  it  is  mentioned  that  the  method  we  used  in  analyzing 
the  eigenvalue  problem  does  not  make  use  of  the  detailed 
structure  of  the  equation  in  the  center  layer  and  is  therefore 
applicable  to  the  magneto-hydrodynamic  case  equally  well. 


•25- 


Acknowledgment 

The  author  wishes  to  thank  Professor  H.  Grad,  Dr,  A. 
Kadish  and  Dr.  G.  Spies  for  several  helpful  discussions. 
Thanks  are  also  due  to  Professor  H.  Grad  for  introducing 
the  author  to  this  subject  matter. 

This  research  was  jointly  supported  by  the  Division 
of  Applied  Mathematics,  Brown  University  and  by  the  U.S. 
A. E.G.  under  Contract  AT( ll-l)-30T7  with  the  Magneto-Fluid 
Dynamics  Division,  Courant  Institute  of  Mathematical 
Sciences,  New  York  University. 


-26- 


References 
*0n  leave  from  the  Division  of  Applied  Mathematics,  Brown 
University 

1.  V.  Aleksin  and  V.  Yashin,  Soviet  Phys.  JETP  13,  787  (1961). 

2.  H.  Grad,  in  Proceedings  of  the  Symposium  on  Electromagnetics 
and  Fluid  Dynamics  of  Gaseous  Plasma  (Polytechnic  Press, 
Brooklyn,  New  York,  1961) 

3.  R.  Kulsrud,  in  Advanced  Plasma  Theory,  ed.  by  M.  N.  Rosenbluth 
(Academic  Press,  1964)  p.  54. 

4.  References  to  some  of  the  earlier  papers  on  this  theory 
can  be  found  in  Refs.  2  and  3. 

5.  K.  Hain  and  R.  Ltist,  Z.  Naturforsch.  13a,  936  (1958). 

6.  H.  Grad,  MHD  Spectrum  and  Low  Shear  Stability,  (to  appear) 

7.  J.  Goedbloed  and  P.  Sakanaka,  A  New  Approach  to  Magneto- 
hydrodynamic  Stability  (to  appear) 

8.  A.  Kadish,  Phys,  Fluids  9,  514  (1966) 

9.  R.  Suydam,  in  Proceedings  of  the  Second  International 
Conference  on  the  Peaceful  Uses  of  Atomic  Energy  (United 
Nations,  Geneva,  1958),  Vol.  31,  p.  157. 


-27- 


Part  2 
LOW  SHEAR  INSTABILITIES  AND  GROWTH  RATES 
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I.   INTRODUCTION 

In  a  recent  work,   a  normal  mode  stability  analysis  has 
been  undertaken  for  a  cylindrical  diffuse  pinch  of  guiding 
center  plasma  to  study  the  kinetic  effect  on  the  stability 
criteria  and  the  growth  rates.   It  is  shown  that  the  problem 
can  be  reduced  to  a  single  second-order  ordinary  differential 
equation  which  may  be  designated  as  the  stability  equation. 
Complete  knowledge  of  the  stability  of  the  pinch  can  only 
come  from  a  thorough  numerical  analysis  on  the  stability 
equation  which  should  also  give  the  growth  rates  associated 
with  any  possible  instabilities.   It  is  difficult  to  get 
general  analytic  results  from  this  stability  equation  because 
of  its  complexity.   For  this  reason  it  seems  worthwhile  to 
go  into  some  special  limiting  situations  where  some  analytic 
results  can  be  drawn.   The  purpose  of  this  work  is  to  present 
some  analytic  results  in  the  low  shear  limit. 

The  results  are  the  line-preserving  instabilities 
considered  in  Sec.  II  for  zero  shear  pinches  and  the 
"semi-local"  modes  in  Sec.  Ill  for  low  shear  pinches.   The 
threshold  and  growth  rates  for  these  instabilities  are 
derived.   It  is  shown  that  for  low  shear  systems  the  large 
end  of  the  "semi-local"  spectrum  joins  the  line-preserving 
spectrum.   This  shows  that  the  line-preserving  instabilities 
have  higher  growth  rate  than  the  "semi-local"  modes  and  this 
justifies  the  zero-shear  approximation  for  a  low  shear  system 
as  far  as  the  maximum  growth  rate  is  concerned.   It  is  also 
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shown  that  when  Suydam-type  instabilities  occur,  the  Suydam 
spectrum  joins  the  small  end  of  the  "semi-local"  spectrum 
which  indicates  at  least  for  the  low  shear  system  the  Suydam 
instabilities  have  lower  growth  rate  and  are  therefore  not 
important.   These  results  are  in  agreement  with  the  corres- 
ponding results  of  the  MHD  theory.    A  comparison  between 
the  present  results  and  the  MHD  results  is  given  in  Sec.  IV, 
where  the  difference  between  the  guiding  center  theory  and 
the  MHD  theory  is  discussed  in  terms  of  the  threshold  and 
growth  rates  of  the  "semi-local"  and  line-preserving 
instabilities . 
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II.   SHEARLESS  DIFFUSE  PINCH 

We  first  consider  the  case  that  the  pinch  is  shearless, 
i.e.  the  rotational  number 

M   =   B^/(rB^)  (1) 

is  a  constant.   It  is  assumed  that  for  the  moment  y  is  equal 
to  a  rational  number  so  that  for  some  suitable  k  and  m  we 
have 

H  =   kB   +  mB,/r   =   0   . 

If  we  use  the  normal  mode  exp(ikz  +  im())  -  iwt)  ,  the 
stability  equation  may  be  written  as 

^  [f^l  *  Gu   =   0  (2) 

where  u  is  the  radial  velocity  component  of  the  plasma  and 
the  coefficients  F  and  G  are  given  by 


rK^(w)NB^[N-E^K^(w)]  (3) 

7      7  7 
po)^  -  O^H  (4) 


G   =   rC^(w,r)  (5) 

w   =   ooB/H  (6) 

2 
and  the  functions  K-,  ,  E,  a   and  C   can  be  found  in  Ref.  1. 

It  suffices  to  note  here  that  if  the  equilibrium  distribution 

functions  are  symmetric  in  the  velocity  v  along  the  field 

lines,  then  the  guiding  center  equations  can  be  shown  to 
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give  rise  to  a  symmetric  operator.    In  such  case,  we  only 
have  to  consider  co  to  be  purely  imaginary  for  the  unstable 
modes.   If,  furthermore,  the  equilibrium  distribution  functions 
are  monotone  decreasing  in  v  ,  then  it  can  be  shown  that 
K, (w)  is  real  and  positive.   Consequently  equation  (2)  is 
non-singular. 

The  boundary  conditions  for  equation  (2)  are  specified 

at  r =  0  and  r=  R,  the  boundary  of  the  plasma.  At  r=  0,  we 
have 

u  =  0     (m^  0)  (7) 

u  =  1     (m=  0)  (8) 

If  the  plasma  is  bounded  by  a  solid  conducting  wall 
at  r=  R,  the  boundary  condition  is  given  by 

u  =  0     at  r  =  R  .  (9) 

And  if  the  plasma  is  separated  from  the  solid  wall  by  a 
vacuum  region,  then  we  have 

a-j^u  +  a^u   =   0    at  r=  R  (10) 

where  a,  and  a^  are  given  in  Ref.  1. 

For  the  case  of  shearless  pinches,  we  set  H  =  0  in  F 
and  G  and  taking  the  limit  w  ->  <»,   The  result  is 

F  =  r[l+  2P2B"2]pY^B2[py2  +  E2(l  +  2P2B"2)]     (n) 

G  =  pY^r  -  IB^B'^  (12) 
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I  =  P^+  P2(1  +  P2B'V  2P^B^r'V2+a2Bjr'^[P2B"2+(2P2  +  P^)(B2  +  2P2)'^]   (13) 

where  y  =  i^  arid  the  dot  denotes  differentiation  with  respect 
to  r. 

It  then  follows  that  for  a  shearless  pinch 

I  >  0  (14) 

is  a  necessary  condition  for  stability.  To  substantiate  this 
claim,  we  only  have  to  consider  small  values  of  y  for  which  F 
and  G  can  further  be  approximated  as 

F   =   rpyVjS^  (15) 

G   =   -Ik^X^  (16) 

2     2     7   2 
where  ^  =  k   +  m  /r  ,  and  let  us  assume  that 

I  <  0 

in  a  subinterval  (r,,r2)  and  I  > 0  in  (o,r^)  and  (r2,R).   We 
write  equation  (2)  as 

d 


^     [pr  ^  37]  -  a2i(B^/B)2u   =   0  (17) 


where  A =  l/y.  Since  X  is  large  when  y  is  small,  (17)  can  be 
solved  by  using  the  WKB  approximation.  Let  us  introduce  the 
following  variables 


r      J^ 

A  |U(t)|^dt  ,  (18) 

^1 


^2      1 
A  [U(t)|^dt  ,  (19) 
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ar)   =   I  k^'(pr)  (20) 

and  first  consider  the  simple  case  where  the  boundary  conditions 
are  given  by  u(o)  =  u(R)  =  0.   In  this  case,  we  require  that 
u(r)  decays  exponentially  as  r  -»■  0  and  r  -^  R  so  that  the 
boundary  conditions  u(o)  =  u(R)  =  0  are  satisfied  with  an 
inaccuracy  which  is  exponentially  small  in  A.   Now,  using 
the  WKB  connection  formulas,   we  have  (1)  to  the  left  of  r^^: 

u   =   a  cos (x^  -  J)  ,  (21) 

(2)  to  the  right  of  T2: 

u   =   a  cos(x_  -  j)    .  (22) 

In  order  that  (21)  and  (22)  agree  with  each  other,  we 
must  have 

A   =   (n  +^)tt/Q  ,  (23) 

where   n   is    an   integer   and 

/2  1 

Q      =      J    lat)|^dt  (24) 

^1       . 


It   then   follows    that    the   growth   rate    is    given  by 

Y      =      Q/[(n4)TT]     ,  (25) 
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which  is  small  if  n  is  taken  to  be  large.   This  shows  that  in 
the  case  where  y  =  B,/(rB^)  is  a  constant  whose  value  is  a 
rational  number,  if  I <  0  for  some  r,  then  there  exists  an 
infinite  number  of  unstable  eigenvalues  going  to  zero. 

If  the  boundary  condition  at  r =  R  takes  the  form  of  (10), 
then  (25)  is  slightly  changed  and  the  result  is 

Y   =   Q/[(n+a)7T]  (26) 

where  the  constant  a  depends  on  the  coefficients  a,  and  a^    in 
(10). 

On  the  other  hand  if  y  =  B  /(rB  )  is  an  irrational  number, 
then  we  can  always  choose  k  and  m  such  that  H  is  small  through- 
out (0,R).   Then  the  foregoing  WKB  analysis  can  still  be 

carried  out  provided  we  restrict  ourself  to  values  of  y  such 

2      2  2 
that  PY   >>  a  H  ,  while  y  itself  is  kept  small.   This 

indicates  that 

I  >  0 

is  also  a  necessary  condition  for  stability  when  \i    is  irrational, 

In  this  case  the  unstable  eigenvalues  are  bounded  away  from 

2      2  2 
the  origin  since  we  have  pY   >>  a  H  . 

We  may  characterize  this  class  of  unstable  modes  by  noting 
that  under  the  condition  H=  0  and  m =  0 ,  the  perturbation  of 
the  magnetic  field  is  aligned  with  the  unperturbed  magnetic 
field  [see  equations  (47)  and  (49)  of  Ref.  1].   For  this 
reason,  these  modes  may  be  termed  the  line-preserving  modes. 
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III.   LOW  SHEAR  STABILITY 

In  the  last  section  the  shearless  systems  are  considered. 
We  shall,  in  this  section,  consider  the  limiting  cases  where 
the  shear  y  is  small,  but  not  necessarily  zero.   We  introduce 
a  small  parameter  e    to  express  the  smallness  o£  the  shear  and 
write 

H   =   mB  (y  +  k/m)   =   eh  (27) 

where  k  and  m  are  specially  chosen  to  make  H  small.   We  again 
consider  small  growth  rate  by  writing 

Y   =   eg  .  (28) 

In  the  limit  e  -^-  0  ,  the  coefficients  F  and  G  in  equation 
(2)  take  the  simple  form 

F   =   e^r[pg2  +  a^H^l/X^  (29) 

G   =   -p(r,YB/H)B^/B2  (30) 

and  p  is  given  by 

p(r,YB/H)  =  P^  +  P2  -Ha^B^/r  +M3(w)+[o2/K^(w)][M2(w)-a2B^]   (31) 

M2(w)  =  B2w232(w)/(rB)  (32) 

MjCw)  =  -w^B32(w)  -  B^[pw^-w^93(w)] (rB^)'^         (33) 

with  the  functions  ^2  ^^^^  ^^'^  ^■zi'^)    given  in  Ref.  1  [equations 

(36)  and  (37)]. 

2 
Since  F  in  (29)  contains  a  small  parameter  e  ,  equation 

(2)  can  again  be  solved  by  the  WKB  approximation  in  a  way 

similar  to  the  analysis  in  Sec.  II.   The  conclusion  is  that 
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a  necessary  condition  for  low  shear  stability  is  that 

p(r,YB/H)   >   0  (34) 

for  all  possible  y   at  all  r,  where  p(r,YB/H)  is  given  by  (31). 
To  substantiate  this,  we  assume  that  p(r,YB/H)  is  negative 
at  some  r  for  Y  =  Y„  =  eg„  and  we  consider  values  of  y  in  a 
neighborhood  of  y    ,    i.e.  y=cg=c{g   +6).   For  6  sufficiently 
small,  we  have 

p(r,gB/h)   <   0 

in  a  subinterval  (r,,r2)  where  r.=r.(6),  i=l,2.   Then 
according  to  the  lowest  order  WKB  approximation,  the  appropriate 
boundary  conditions  at  r=0  and  r=R  can  be  satisfied  if 

e'^L(6)   =   nTT  +  a  (35) 

for  n =  0  ,±1  ,±2  , . . .  where  a  is  the  same  constant  used  in  (26) 
and  L(6)  is  given  by 


^2  - 

LC6)  =  j  |p(r,gB/h)k2/[(a2h2+pg2)r]|2dr     (36) 

^1 

where  it  is  recalled  that  g  =  g  +5  and  r,  =  r,  (6)  ,  T2=   r^  (6)  . 

For  sufficiently  small  6,  then  (35)  can  be  written  as 
6L'  (0)  =  e(})(n,e)  (37) 

())(n,e)  =  mr  +  a  -  -L(0)  (38) 
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where  L'(0)  is  the  derivative  of  L(6)  with  respect  to  6  at 

es 

(39) 
We  write  (37)  as 


6  =  e[(l)(n^,£)  +  Z-n]  (40) 

which  determines  6  as  1  =  0 ,±1 ,±2 , .  .  .  .   It  is  noted  that  ei 

must  be  small  so  as  to  keep  6  small  as  it  is  assumed. 

This  shows  that  if  for  some  y  >  0  we  have  p(r,Y  B/H)  <  0 


at  some  r,  then  it  is  always  possible  to  choose  a  y  which 
differs  : 
;iven  by- 


differs  from  y   in  an  amount  of  the  order  e      and  which  is 


with  6  given  by  (40)  such  that  the  boundary  conditions  at 
r=0  and  r=R  are  satisfied.   Consequently,  the  system  is 
unstable.   Furthermore,  it  follows  from  (40)  that  the  spacing 

of  the  unstable  eigenvalues  is  of  order  e  .   These  instabilities 

2 
correspond  to  the  "semi-local"  modes  in  the  MHD  diffuse  pinch. 

For  these  instabilities  the  maximum  growth  rate  is  given  by 

y     =   Max    {y:  p(r,yB/H)  <  0}  (41) 

^^"^   0<r<R 

with  an  error  of  e  -order. 

It  is  worthwhile  to  note  that  as  yB/H  ^  °° , 

lim  p(r,yB/H)   =   I  (42) 
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where  I  is  given  by  expression  (13).   This  seems  to  suggest 
that  i£  I  <  0  at  some  r  then  expression  (40)  would  imply  that 
the  maximum  growth  rate  is  unbounded.   What  this  actually  means 
is  that  the  growth  rate  y  is  no  longer  of  order  e  when  I  becomes 
negative.   To  see  this,  we  return  to  (37)  and  consider  e  fixed 
and  let  g^  =  y J z   become  large.   Then  the  coefficient  L' (0)  of 
6  becomes  increasingly  small  since 

L'(0)  -^   1/go  , 
which  shows  that  in  order  for  6  to  remain  small  g  must  remain 
much  smaller  than  1/e.   In  fact,  for  large  y  (i.e. 
PY^  >>  a^H^),  we  can  neglect  a^H^.   This  will  lead  us  back  to 
the  result  of  Sec.  II  and  the  resultant  growth  rate  is  then 
given  by  (25).   In  other  words,  when  I  becomes  negative,  the 
large  end  of  the  closely  spaced  spectrum  of  the  "semi-local" 
modes  joins  the  spectrum  of  the  line-preserving  modes  examined 
in  Sec.  II. 

On  the  other  hand  as  yB/H  -+  0 ,  we  have 

p(r,o)  =  P^+  P2+  a2(B^/r)[l  -  0^X^(0)]  . 

Now  if  p(r,o)  <  0  at  some  r,  then  there  exist  unstable 
semi-local  modes  whose  eigenvalues  are  given  by  (40)  as 


where  6  can  be  calculated  from  (37)  and  must  remain  positive. 
This  shows  that  y  is  bounded  away  from  the  origin  by  a  distance 
of  order  e  .   Now  if  the  guiding  center  Suydam  condition   is 
violated;   i.e.  if  H=  0  at  r=  r   and  if 
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S  =  p(r,o)  +  ^^^liu/v)^    <   0 


at  r=r  ,  then  there  appear  infinitely  many  Suydam  modes 
whose  eigenvalues  go  to  zero  with  the  origin  as  an  accumula- 
tion point.   This  shows  that  when  the  low  shear  system  is 
Suydam  unstable  the  small  end  of  the  semi-local  spectrum 
joins  the  Suydam  spectrum. 
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IV.   CONCLUDING  REMARKS 

In  the  preceding  sections  we  have  studied  the 
line-preserving  instabilities  in  a  zero-shear  diffuse  pinch 
and  the  semi-local  instabilities  in  a  low  shear  diffuse  pinch. 
It  is  shown  that  the  large  end  of  the  semi-local  spectrum 
joins  the  line-preserving  spectrum  and  the  small  end  joins 
the  Suydam  spectrum.   This  is  in  agreement  with  the  MHD 
result  of  Grad.    This  shows  that  the  zero-shear  approximation 
is  a  good  model  for  low  shear  systems  in  finding  out  the 
maximum  growth  rate  given  by  the  line-preserving  unstable 
modes . 

This  suggests  that  a  possible  estimate  for  the  maximum 
growth  rate  for  a  low  shear  diffuse  pinch  may  be  obtained 
from  (26)  as 

YMax   '^  ?  j  Ik'l/(pr)l^dr  (45] 

^1 
Since  I,  given  in  (13),  contains  only  P-|^,  ^2    ^^^  ^* 

the  threshold  and  the  growth  rate  of  the  line-preserving 

instability  do  not  depend  on  the  details  of  the  distribution 

functions.   In  the  scalar  pressure  case  we  have  from  (13) 

I  =  P(2 +  P/B^)  + 6(B^/r) (P/B^) (B^+P)/(B^+2P)  (46) 

while  the  corresponding  MHD  theory  gives 

I  =  2P  +  4(B^/r)pa^/(B2  +  pa^)  (47) 

where  a   =  9P/3p;  i.e.  a  is  the  sound  speed.   For  low-3 
plasma  (46)  and  (47)  can  be  made  to  agree  by  taking  the  ratio 
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o£  specific  heats  C  /C^  =  3/2.   But  in  high-3  plasma  no 
such  simple  relation  exists  between  (46)  and  (47). 

In  the  cases  of  scalar  pressure  where  B,  is  much  smaller 
than  the  longitudinal  field  component  B  ,  the  Suydam 
instability  condition  in  (44)  implies  the  line-preserving 
instability  I  <  0  where  I  is  given  by  (46).   This  reinforces 
the  importance  of  the  line-preserving  instability  over  the 
Suydam  instability  for  low  shear  system.   In  other  words, 
if  a  low  shear  system  with  B   <<  B   is  Suydam  unstable,  then 
it  has  line-preserving  unstable  modes  which  give  the  maximum 
growth  rate. 
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